Geometric strong segregation theory for compositionally asymmetric diblock 

copolymer melts 
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We have identified the effect of the Wigner-Seitz ceil geometry in the strong segregation limit of 
diblock copolymer melts with strong composition asymmetry. A variational problem is proposed 
describing the distortions of the chain paths due to the geometric constraints imposed by the cell 
shape. We computed the geometric excess energies for cylindrical phases arranged into hexagonal, 
square, and triangular lattices and explicitly demonstrated that the hexagonal lattice has the lowest 
energy for a fixed cell area. 
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Block copolymers are a well-known class of smart 
materials that can produce a wide variety of com- 
plex equilibrium microstructuresii^i^i^. Since the late 
70's, these systems have received significant atten- 
tion by theorists, and many of the copolymer system 
phases are now well explained on the basis of energy 
minimization arguments^'^'^ i^'^'^^i^^'^^i"'^^ . Further ad- 
vances in computational techniques allowed studies of 
the phase behavior in block copolymer systems with 
the help of direct numerical solution of models which 
explicitly incorporate statistical mechanics of polymer 
chains, providing a connection between the observed mi- 
crostructures and the underlying microscopic material 
parametersi^ii^iiiii^ii^i2ai2ii2^. More recently, a re- 
newed interest in block copolymer systems was stimu- 
lated by the studies of geometrically constrained systems, 
such as block copolymers confined to the surface of a sub- 
strate or the inside of nanopores^^i^li^. Numerical stud- 
ies of these systems showed a rich variety of microstruc- 
tures with very intricate geometrie a^^i^^'^^i^^ . 

An important regime in which the self-organizing be- 
havior of block copolymer systems becomes especially 
pronounced is the strong segregation limit, in which 
monomers of different types segregate almost completely 
into non-overlapping regions of space. In the case of di- 
block copolymers, the first theory of the strong segre- 
gation limit was proposed by Semenov, who computed 
the phase coexistence boundaries for several types of mi- 
crophases with simple geometrieai^. His results were es- 
sentially corroborated by the computational studies of 
self-consistent mean-field models and extended under a 
number of simplifying assumptions to include more exotic 
phasesi^ii^iiiiii^. Nevertheless, the original Semenov's 
theory crucially relies on approximating the Wigner-Seitz 
cell of the corresponding periodic structure by a disk of 
the same area in the case of cylindrical phases, or a ball 
of the same volume in the case of spherical phases. Thus, 
Semenov's theory neglects the effect of the cell geometry 
and, therefore, cannot distinguish between the structures 
which are characterized by different unit cell types (as, 



e.g., body-centered cubic versus face-centered cubic lat- 
tices of spheres) . This point becomes even more critical 
when considering the effects of confinement, since the do- 
main shape must be truly important in such problems. 

In this paper, we propose an extension of the strong 
segregation theory for block copolymer systems which 
includes the effect of geometric constraints on the chain 
configurations. Specifically, we investigate the case of 
diblock copolymers with strong composition asymmetry, 
in which most of the "unpleasant" features of the strong 
segregation limif^ are under control, allowing us to con- 
centrate on essentially the only remaining issue of the 
effect of the geometry. We formulate a variational prob- 
lem that gives the excess energy due to geometric factors 
and, in particular, allows to discriminate between differ- 
ent lattice types with the same unit cell volume in the 
case of periodic microstructures. To illustrate the lat- 
ter point, we explicitly compute the excess energies for 
cylindrical phases on three fundamental two-dimensional 
lattices and demonstrate that, as intuitively and experi- 
mentally expected, the hexagonal lattice has the lowest 
energy. 

Consider a system of linear polymer chains consist- 
ing of N monomers of type A bonded covalently with 
JN monomers of type B, with / <C 1. If each A- and 
B-monomer has excluded volume v, then / is basically 
the volume fraction of the B-monomer. We introduce 
the Flory interaction parameter the Kuhn statistical 
length h, and the root-mean-square end-to-end distance 
R ~ b^/N. The confinement domain (or the Wigner-Seitz 
cell for periodic structures) is denoted by C K"^. 

In the strong segregation regime xN ^ 1 the A- and B- 
monomers locally segregate into disjoint subsets Ha and 
Hb of O, with a sharply defined interface F = dftA<^dftB 
containing the A-B jmictions. Based on this observation 
and the Gaussian chain model, Semenov computed the 
free energy of the system as a sum of three contributions: 



interface 



(1) 
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Here, ^interface is the interfacial energy given by 



-Fintorfaco ^ j 



(2) 



where (Tq is a dimensionless parameter of order 1, ks 
is the Bohzmann constant and T is temperature. Next, 
-fcore is the energy of the smah B-monomer core, assumed 
to be radially symmetrioi^ii^iii^: 



/ z^(r)dr, 



(3) 



where z(r) is the distance from r to F. Finally, Fcorona 
is the energy of the corona composed of the A-monomers 
filling flA- Both calculations of -Fcorc and -Fcoiona rely on 
the strongly stretched Gaussian chain assumption. How- 
ever, the computations differ crucially, since the surface 
F, as seen from Ha, is convex, while from the side of f2s it 
is concave. A parabolic brush assumption can be used to 
compute i^coro, while one needs to take into account the 
exclusion zone to compute i^^.o^onET^'^'^'^'^'^'^- In 
the original Semenov's theory, the domain is replaced 
with a cylindrical or spherical region of the same volume, 
and the Alexander-de Gennes brush assumption^ii^ is 
used to compute the energy. This, however, does not af- 
fect the leading order contribution to the energy at small 
volume fractions, since the energy is dominated by the 
singularity near the interface^Si^iHi^. 

We will now compute the corona energy, taking also 
the geometric effects into account. We start with the 
self-consistent mean-field theory in which each polymer 
molecule is treated as a Gaussian chain with the position 
To S F of the A-B junction uniformly distributed on the 
interface (the latter is justified for / ^ 1). Then the 
corona free energy is (up to an additive constant) 



ksT 



Here, ri, . . . , rjv denote the positions of the A-monomers, 
ro is the position of the A-B junction, v is the density of 
the A-B junctions on F, and Ha is the Gaussian chain 
Hamiltonian: 



ikeT 
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^(r„ -r„_i)^ 
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where ipA is the self-consistent field (a Lagrange multi- 
plier) enforcing the average monomer density to be 
everywhere in ^Ia- 

In the strong segregation limit the chains are highly 
stretched, i.e. we have Ir^r — ro| ^ i?g, so one would 
naturally want to use the method of steepest descent to 
evaluate the integral in However, to proceed further 



we note a general difficulty that in the strong segregation 
limit the integral in ^ is not dominated by the global 
minimizer of Ha, the fluctuations of the end-point po- 
sitions ym actively contribute to the free energy of the 
chain&iS. Nevertheless, when F is convex, as seen from 
VLa , an exclusion zone must form around b which is free 
of the chain ends^i^. Moreover, for / ^ 1 this exclusion 
zone must occupy most of Q.a, pushing chain ends close 
to the cell boundary. In the case when f2 is a disk an ex- 
act solution to the problem shows that the exclusion layer 
extends to the fraction of 2/7r ~ 0.64 of the disk radius, 
and in fact the majority of the chain ends are located 
within about 6% of the outer boundary^. For a sphere 
the distribution of chain ends is even tighter, with the 
dead layer extending to about 0.76 of the radius, with 
the majority of the ends within about 4% of the outer 
boundary^'^. Therefore, for sufficiently small values of / 
a very good approximation to the problem with an ex- 
clusion zone should be given by the Alexander-de Gennes 
brush^iiS^, in which all chain ends are assumed to lie on 
the outer boundary 90^. In the following, we adopt this 
approximation to eliminate the need to deal with the pre- 
cise chain end statistics. We also note that this assump- 
tion is expected to be asymptotically exact when the free 
end of the A-chains is capped by sticky end-groups^i^ 
or by a short block of C-monomers which is immiscible 
with either A- or B-blocks. 

Under the assumption of Alexander-de Gennes brush 
for the corona, the integral in Q is dominated by the 
minimizers of Ha with Tq G F fixed and r^v restricted to 



-^corona — ^ 

Now, introducing 
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and then passing to continuous chains: r„ ~ r{n/N), 
where r : [0, 1] SIa are continuous paths, we can write 
the corona energy as 



iksTv 



r Jo 



dt dA, 



where each path r{t) minimizes 

2 



S = 



Uir) }dt 



(8) 



(9) 



with fixed endpoints. 

Note that to be in the mechanical equilibrium the 
chains must come out normally from the interface F. 
Then, from the constant monomer density requirement 
near F one can get the initial conditions for the minimiz- 
ers r(i), given a potential U{r) enforcing the constraint: 



r(0) = ro. 



dr(0) 
dt 



Nvvn{ro), 



(10) 



3 




where n is the outward (from fls) normal to F at rg G T. 

We now point out a mechanical analogy, according to 
which r{t) can be interpreted as the trajectory of a point 
particle with unit mass in M.^ moving under the action 
of potential energy U. The function S plays the role of 
the action^. Therefore, the equation of motion for r{t) 
becomes simply 



-VC/(r). 



(11) 



Note that the initial condition in ([TO]) then uniquely 
determines the point at which the trajectory r{t) hits 
dfl. We can also easily write down the corresponding 
Hamilton- Jacobi equation'^^: 



(12) 



where Ug is a constant which a posteriori turns out to be 
independent of the initial point Tq of the trajectory that 
passes through r. 

Note that we still need to determine the self-consistent 
field, now given by U, enforcing the constant monomer 
density constraint, which in terms of r{t) becomes 



vN 



JT Jo 



v{t))dtdA = V' 



(13) 



where 5{r) is the three-dimensional delta-function, and 
we assumed that the family of trajectories v{t) foliates 
VIa- On the other hand, observe that it is, in fact, suf- 
ficient to find only the action S appearing in p2|) . The 
Hamiltonian structure of equation of motion ITTI) imposes 
certain restrictions on the possible trajectories r(t), in 
particular, it forces the dynamics of r to be a gradient 
flow. In fact, it is easy to see that this gradient flow 
also has to be divergence-free. Indeed, consider a tube 
formed by trajectories originating on some closed curve 
in r enclosing an area A, and write down the total num- 
ber M of A-monomers contained in the cylinder between 
t = to and t = to + T cross-sections of that tube. One 
easily gets M = vNAt, which is clearly independent of 
tf). Hence, differentiating this quantity with respect to 



to, we see that the total flow in/out of the cylinder along 
the trajectories must equal zero. In view of arbitrariness 
of to,T and A, we must have V • (dr/dt) = in ^Ia- 

The arguments above immediately imply that the ac- 
tion S must be a harmonic function: 



AS* = in nA, 



n • VS" = Nvi^ on T, 



(14) 



where A is the Laplacian, and we also used pOj) . On the 
other hand, the boundary data on dfl must be chosen in 
an unusual way: every trajectory starting on T at t — 
and flowing up the gradient of S must reach dQ at t = 1. 
This condition can also be reformulated as: 



1 on every field line of S. 



|V5(r)| 

It is also easy to see from (|T3)) that 





dv 




u: 




dtdA= [ 




'dt 


Jn 



|V5pdr. 



(15) 



(16) 



With this, the expression for the corona energy becomes 
SksT 



D[Si 



D[S] = [ \VS\^dr. (17) 



Let us note that existence and uniqueness of solutions 
to the proposed problem is not guaranteed a priori for 
any given Q and F. In view of (|17p. one should, in fact, 
be interested in minimizers of D[S] which also satisfy 
the condition in (|15p . Notice that the location of F rel- 
ative to Q is also part of the minimization problem. A 
radial solution (presumably, the unique minimizer) triv- 
ially exists when Qb smd are concentric balls. One 
would then expect from perturbative considerations that 
this solution should persist when dfl is slightly distorted 
away from a sphere. 

In view of the assumption / ^ 1, the solution of (fT4|) 
coincides to the leading order with that of 



AS= \n\6{r). 



(18) 



Indeed, if 0^ is a ball of radius R centered at the origin, 
then AttR^vNv = \Q\ to the leading order in /, where 



4 



denotes the volume of 17 (same argument applies in two 
dimensions). From this, one can see that since S should 
behave as the free space Green's function of the Laplacian 
near the origin, D[S] will diverge as i? ^ 0. The leading 
order singular term will only depend on |n| and not the 
shape of fl and is precisely what was calculated by Se- 
menov for the corona energy On the other hand, for 
small but finite values of R the solution will also contain 
an excess energy associated with the geometry of fl. 

We applied our variational procedure to compute the 
excess energy in the case of two-dimensional hexagonal, 
square, and triangular lattices of straight cylinders Qb 
of radius R and unit height. We note first that the exact 
solution of the problem in a coaxial cylinder of the same 
volume |r2| gives straightforwardly 

D,= ^^lnr\ (19) 

To compute the excess energy D — Dq for Wigner-Seitz 
cells corresponding to the considered lattices, we imple- 
mented a finite element-based minimization algorithm to 
find minimizers of D satisfying (fT5|)i2.. The minimizing 
trajectories in cells whose area is normalized to unity are 
presented in Fig. [TJ From dimensional arguments, we 
find that for / ^ 1 we have 

D^Do^Cn\n\^, (20) 

where the dimensionless constant Ca depends on the ge- 
ometry of the cell only. Numerically, we found Chcx — 
0.00922, Csq ~ 0.0453, and Ctriang ^ 0.179 for the hexag- 
onal, square, and triangular lattices, respectively. Note 
that for fixed values of R and / both the interfacial en- 
ergy Fintorfaco, thc corc energy i^coro, and the leading- 
order corona energy -Fcorona obtained from ([T7]) with D 
replaced by Dq are the same. Therefore, to compare the 
energies of different geometric arrangements of the B- 
domains, one needs to compare the excess energies. From 
our calculation above we can immediately conclude that 
among the considered types of lattices of cylinders with 
the same radius R and volume fraction / the hexago- 
nal lattice is the most energetically favorable in the limit 



/ ^ 0, an intuitively expected result which is put on a 
rigorous footing by our computations. Let us point out 
that our approach should also be applicable to spheri- 
cal phases to help identify the minimizer among differ- 
ent types of three-dimensional lattices. We note that 
the answer to this question in the strong segregation 
limit lies beyond the scope of Semenov's theory^^ and its 
extensionsi^i^. Let us also point out that the method 
of Refsii^ii^ cannot be applied here, since it ignores the 
effect of the exclusion zone. 

Let us note that our calculation is akin to the one 
performed by Fredrickson^, who estimated the excess 
energy due to geometric factors for a hexagonal Wigner- 
Seitz cell in the strong segregation limit. Fredrickson 
used linear elasticity and the Alexander-de Gennes as- 
sumption to study the extra contribution to the elastic 
energy of the corona due to chain distortions. His re- 
sult, however, differs from ours quantitatively. In partic- 
ular, we find the excess energy obtained by us is greater 
than thc one obtained by Fredrickson by a factor of 1.5. 
We attribute this discrepancy to strong chain distor- 
tions, which invalidate the linear elasticity approxima- 
tion. Thus, at small / the excess energy due to geometry 
of the Wigner-Seitz cell may have a larger contribution 
than previously expected. 

To conclude, we have developed a variational charac- 
terization of the leading geometric corrections to the Se- 
menov's strong segregation theory in the case of strong 
composition asymmetries. Our theory thus should be 
able to account for the effect of the confinement geome- 
try on microstructures consisting of small droplets of the 
minority species and, in particular, help identify the equi- 
librium lattice configurations of these droplets, as was 
explicitly demonstrated in the case of cylindrical phases. 
Perhaps more importantly, our theory provides a new 
way to study questions of metastability and instability 
of nonequilibrium copolymer microstructures under ex- 
ternal perturbationsiii^i^i^. 
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